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ABSTRACT 

We  analyze  a  sample  of  58  multiwavelength,  Very  Long  Baseline  Array  observations  of  active  galactic  nuclei 
(AGNs)  to  determine  their  scattering  properties.  Approximately  75%  of  the  sample  consists  of  AGNs  that  exhibit 
centimeter-wavelength  intraday  variability  (interstellar  scintillation),  while  the  other  25%  do  not  show  intraday  vari¬ 
ability.  We  find  that  interstellar  scahering  is  measurable  for  most  of  these  AGNs,  and  the  typical  broadening  diameter  is 
2  mas  at  1  GHz.  We  find  that  the  scintillating  AGNs  are  typically  at  lower  Galactic  latitudes  than  the  nonscintillating 
AGNs,  consistent  with  the  scenario  that  intraday  variability  is  a  propagation  effect  from  the  Galactic  interstellar  me¬ 
dium.  The  magnitude  of  the  inferred  interstellar  broadening  measured  toward  the  scintillating  AGNs,  when  scaled  to 
higher  frequencies,  is  comparable  to  the  diameters  inferred  from  analyses  of  the  light  curves  for  the  more  well-known 
intraday  variable  sources.  However,  we  find  no  difference  in  the  amount  of  scattering  measured  toward  the  scintillating 
versus  nonscintillating  AGNs.  A  consistent  picture  is  one  in  which  the  scintillation  results  from  localized  regions 
(“clumps”)  distributed  throughout  the  Galactic  disk,  but  that  individually  make  little  contribution  to  the  angular  broad¬ 
ening.  Of  the  58  AGNs  observed,  37  (64%)  have  measured  redshifts.  At  best,  a  marginal  trend  is  found  for  scintillating 
(nonscintillating)  AGNs  to  have  smaller  (larger)  angular  diameters  at  higher  redshifts.  We  also  use  our  observations  to 
try  to  constrain  the  possibility  of  intergalactic  scattering.  While  broadly  consistent  with  the  scenario  of  a  highly  turbulent 
intergalactic  medium,  our  observations  do  not  place  significant  constraints  on  its  properties. 

Subject  headings:  galaxies:  active  —  galaxies:  ISM  —  galaxies:  jets  —  galaxies:  nuclei  —  ISM:  structure  — 
quasars:  general  —  radio  continuum:  galaxies  —  surveys 


1.  INTRODUCTION 

There  is  now  compelling  evidence  that  the  intraday  variability 
(IDV)  phenomenon — intensity  variations  on  hour  timescales  at 
centimeter  wavelengths  in  compact,  flat-spectrum,  extragalactic 
sources  (e.g.,  Heeschen  1984;  Quirrenbach  et  al.  1992) — is  of 
extrinsic  origin.  Density  fluctuations  in  the  interstellar  medium 
(ISM)  induce  refractive  index  fluctuations  that,  when  combined 
with  the  relative  motions  of  the  scattering  medium  and  the  Earth, 
produce  intensity  variations  or  scintillations.  The  evidence  for 
this  conclusion  is  both  differences  in  the  variability  pattern  arrival 
times  at  widely  spaced  radio  telescopes  and  annual  cycles  in  the 
variability  characteristics  for  various  IDV  sources  (Jauncey  et  al. 
2000, 2003,  2006;  Jauncey  &  Macquart  2001;  Rickett  et  al.  2001; 
Dennett-Thorpe  &  de  Bruyn  2002,  2003;  Bignall  et  al.  2006). 

In  order  to  exhibit  interstellar  scintillations  (ISS),  a  source  must 
contain  a  sufficiently  compact  component  (analogous  to  “Stars 
twinkle,  planets  do  not”)  such  that  its  angular  diameter  is  com¬ 
parable  with  or  smaller  than  the  size  of  the  first  Fresnel  zone  of 
the  scattering  screen,  i.e.,  of  order  tens  of  microarcseconds  at  fre¬ 


quencies  near  a  few  gigahertz  (e.g..  Walker  1998).  The  absence 
of  ISS  in  a  source  could  be  either  an  intrinsic  or  extrinsic  effect. 
Active  galactic  nuclei  (AGNs)  might  have  a  range  of  intrinsic  di¬ 
ameters,  in  which  case  only  the  most  compact  would  exhibit  ISS. 
Alternately,  interstellar  density  fluctuations  produce  a  rich  range 
of  observable  phenomena  (Rickett  1990),  of  which  scintillations 
are  only  one  manifestation.  Angular  broadening  along  the  line 
of  sight,  due  to  either  multiple  ionized  media  or  an  extended 
medium,  could  produce  apparent  diameters  of  AGNs  sufficiently 
large  that  the  AGNs  would  not  display  ISS. 

Consistent  with  this  requirement  for  a  compact  (^10  /ras)  com¬ 
ponent,  Ojha  et  al.  (2004)  have  compared  AGNs  that  display  ISS 
with  those  that  do  not,  and  find  that  the  scintillating  AGNs  typi¬ 
cally  are  more  core-dominated  on  a  milliarcsecond  scale  than  the 
nonscintillating  AGNs.  This  result  is  striking  given  that  their  ob¬ 
servations  compared  the  source  structure  on  milliarcsecond,  not 
microarcsecond,  scales.  However,  the  observations  of  Ojha  et  al. 
(2004)  were  at  the  single  frequency  of  8.4  GHz.  This  frequency 
is  sufficiently  high  that  interstellar  scattering  effects  on  most  lines 
of  sight  through  the  ISM  would  not  be  detectable  on  terrestrial 
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Notes. — Sources  are  indicated  to  be  either  scintillators  (Y),  nonscintillators  (N),  or  sources  observed  to  vary  once  and  presumed  to  be  scintillators  (Yj);  a:  is  the 
spectral  index  for  the  frequency  dependence  of  the  intrinsic  diameter,  eq.  (1);  9^  and  Oi  are  the  inferred  scattering  and  intrinsic  diameters,  respectively;  N  and  ^re  the 
number  of  data  and  chi-square  in  the  fit  for  eq.  (1);  and  ^ne200i  is  the  predicted  interstellar  scattering  diameter  from  the  NE2001  model. 

^  Both  and  dependences  for  the  intrinsic  diameter  produced  equal  x^- 
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baselines,  nor  could  they  have  exploited  the  wavelength  depen¬ 
dence  for  scattering  to  separate  its  effects  from  the  intrinsic  diam¬ 
eters  of  the  sources. 

The  AGNs  observed  by  Ojha  et  al.  (2004)  were  drawn  from 
the  Micro-Arcsecond  Scintillation-Induced  Variability  (MASIV) 
survey  (Lovell  et  al.  2003).  The  MASIV  survey  is  a  large  variabil¬ 
ity  siuwey  of  the  northern  sky  with  the  primary  goal  being  the  con¬ 
struction  of  a  large  sample  of  scintillating  AGNs.  The  survey  used 
the  Very  Large  Array  ( VLA)  at  5  GHz  in  a  multiarray  mode  and 
has  yielded  scintillation  information  on  over  500  AGNs,  of  which 
over  half  have  been  found  to  be  scintillating  (Lovell  et  al.  2007). 

Ojha  et  al.  (2006)  presented  multifrequency  observations  of 
a  subset  of  MASIV  sources.  Their  observations  were  designed 
to  create  a  sample  of  sufficient  size  to  compare  and  contrast  the 
scattering  behavior  of  scintillating  and  nonscintillating  AGNs. 
This  paper  reports  the  analysis  of  those  observations. 

2.  SOURCE  SAMPLE 

Our  sample  consists  of  58  AGNs,  observed  in  two  different  pro¬ 
grams.  The  first  subset  consists  of  49  AGNs  from  the  MASIV 
survey,  observed  over  3  days  with  the  Very  Long  Baseline  Array 
(VLBA)  in  2003  February.  At  the  time  of  the  observations, 
approximately  half  of  these  49  AGNs  were  classified  as  highly 
variable  MASIV  sources,  with  scintillation  indices  larger  than 
2%;  the  other  half  were  classified  as  nonscintillators,  with  no 
scintillation  index  larger  than  0.2%.  These  sources  were  chosen 
without  regard  to  their  Galactic  latitude  or  longitude.  Since  our 
observations,  however,  further  MASIV  observations  and  anal¬ 
ysis  shows  that  many  of  the  AGNs  identified  originally  as  non¬ 
scintillating  are  in  fact  scintillating  ( Lovell  et  al.  2007 ).  The  number 
of  recognized  scintillating  AGNs  is  now  46,  and  the  number 
of  nonscintillating  sources  is  now  1 3  (ratio  of  3 : 1 ).  The  subset  of 
AGNs  from  the  MASIV  survey  were  observed  at  0.33,  0.61,  and 
1 .6  GHz,  with  additional  observations  at  2.3  and  8.4  GHz.  In  total, 
1 94  VLBA  images  of  the  49  MASIV  extragalactic  radio  sources 
at  up  to  seven  observing  frequencies  were  obtained  (Ojha  et  al. 
2006 ).  Additional  data  were  obtained  from  the  United  States  Naval 
Observatory  (USNO)  Radio  Reference  Frame  Image  Database 
(RRFID)^  and  the  literature. 

The  observations  of  these  49  AGNs  were  acquired  by  cycling 
through  the  sources  so  as  to  increase  the  u-v  plane  coverage.  Typ¬ 
ical  times  on  source  range  from  10  minutes  at  the  higher  frequen¬ 
cies  to  25  minutes  at  the  lower  frequencies.  Typical  noise  levels 
were  within  a  factor  of  2-3  of  the  expected  thermal  noise  limits. 

Ojha  et  al.  (2006)  fit  Gaussian  component  models  to  the  vis¬ 
ibility  data  of  the  sources,  using  the  images  as  guides.  If  more  than 
one  component  was  required  to  model  a  source  at  a  given  fre¬ 
quency,  the  most  compact  component  was  identified  as  the  “core,” 
as  the  most  compact  component  will  be  the  one  for  which  scat¬ 
tering  effects  will  be  most  apparent.  The  most  compact  compo¬ 
nent  is  frequently  the  brightest  one.  For  the  few  sources  where 
this  does  not  hold  strictly  at  all  frequencies,  the  compact  and 
bright  component  that  could  consistently  be  identified  as  the  same 
at  all  frequencies  was  identified  as  the  core,  e.g.,  J07 13+4349, 
where  the  northernmost  component  is  identified  as  the  core  even 
though  it  is  not  the  brightest  component  at  8.4  GHz  (Ojha  et  al. 
2006). 

The  second  subset  consists  of  nine  AGNs  chosen  from  those 
used  to  define  the  International  Celestial  Reference  Frame  (Johnston 
etal.  1995;Maetal.  1998).  The  initial  motivation  was  to  use  these 
AGNs  to  search  for  scattering  resulting  from  the  interstellar  media 

’  See  http://rorf.usno.navy.mil/rrfid.shtml. 


Fig.  1 . — Two  examples  of  fitting  results.  For  both  plots,  circles  show  measured 
diameters,  arrows  indicate  upper  limits,  the  solid  line  indicates  the  fit  of  eq.  (1)  to 
the  observations,  the  dashed  line  indicates  the  inferred  scattering  diameter,  and  the 
dotted  line  indicates  the  inferred  intrinsic  diameter.  Uncertainties  on  the  angular 
diameters  are  plotted,  but  in  many  cases  are  comparable  to  the  size  of  the  symbol. 
Top:  The  source  J2022+6136,  for  which  a  relatively  large  scattering  diameter, 
3.1  mas  at  1  GHz,  is  inferred.  Bottom:  The  source  J0745+1011,  for  which  a  rela¬ 
tively  small  scattering  diameter,  0.6  mas  at  1  GHz,  is  inferred. 

of  galaxies  along  their  lines  of  sight.  The  sources  were  chosen  to 
be  (1)  at  Galactic  latitudes  \b\>  45°;  (2)  strong,  with  cm  ~  1  Jy; 
and  (3)  compact  and  dominated  by  a  single  component  (Fey  & 
Chariot  1997).  Observations  of  these  nine  sources  were  conducted 
with  the  VLBA  on  2001  Febraary  17  and  18at0.33  and0.61  GHz. 
Calibration,  imaging,  and  model  extraction  were  performed  in  a 
marmer  similar  to  that  used  by  Ojha  et  al.  (2006). 

Because  both  subsets  involved  observations  at  0.33  and 
0.61  GHz,  the  observations  were  generally  carried  out  at  night  so 
that  the  sources  had  large  solar  elongations.  Indeed,  observing  at 
large  solar  elongation  was  an  explicit  criterion  in  scheduling  the 
observations  for  the  second  subset.  As  a  result,  the  smallest  elon¬ 
gation  for  any  source  is  75°,  and  the  typical  elongation  is  approx¬ 
imately  130°. 

For  the  present  analysis,  we  have  used  the  core  component 
models  from  these  two  observing  programs,  augmented  by  mea¬ 
surements  from  the  literature.  All  sources  have  angular  diameters 
measured  at  at  least  three  frequencies,  and  some  sources  have 
measured  angular  diameters  at  as  many  as  seven  frequencies.  See 
Table  2  of  Ojha  et  al.  (2006).  While  the  selection  criteria  for  the 
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Fig.  2. — Distribution  of  sources  observed  as  a  function  of  Galactic  coordinates.  The  Galactic  anticenter  is  at  the  center  of  the  plot,  and  longitude  increases  to  the  left. 
Stars  indicate  sources  that  scintillate;  squares  indicate  nonscintillating  sources.  The  size  of  the  symbol  is  qualitatively  proportional  to  the  inferred  scattering  diameter. 


two  subsets  differed,  the  sources  were  treated  identically  in  the 
following  analysis. 

3.  ANALYSIS 

Angular  broadening  is  manifested  as  an  observed  angular  di¬ 
ameter  that  scales  approximately  as  1^.  We  fit  the  measured  an¬ 
gular  diameters  to  the  functional  form 

0^  =  +  (Oyf,  (1) 

where  0s  and  0,  are  the  scattering  and  intrinsic  (FWHM)  diam¬ 
eters  of  the  AGN,  respectively,  at  the  fiducial  frequency  of  1  GHz. 
We  found  fhe  best-fitting  values  for  0s  and  0i  in  a  minimum  x  ^ 
sense.  We  considered  both  x  =  0  (i.e.,  frequency-independent 
intrinsic  diameter,  for  a  flat-spectrum  source)  and  x  =  —  1  (i.e., 
frequency  scaling  for  a  single  incoherent  synchrotron  component) 
and  selected  the  value  of  x  that  produced  the  lower 

As  a  motivation  for  the  use  of  equation  (1),  as  well  as  antic¬ 
ipating  later  discussion,  we  begin  by  considering  a  crude  approx¬ 
imation  to  equation  (1).  For  the  sources  having  measured  angular 
diameters  at  both  0.33  and  1.6  GHz,  we  have  assumed  a  simple 
power-law  scaling  for  the  angular  diameter,  0  ex  and  solved 
for  /?.  We  chose  0.33  GHz  because  the  frequency  dependence  of 
scattering  means  that  it  will  be  the  most  sensitive  to  scattering;  we 
chose  1 .6  GHz  as  the  second  frequency  as  an  attempt  to  balance 
between  having  a  sufficiently  large  frequency  dynamic  range  so  as 
to  obtain  a  robust  estimate  of  (3,  but  not  having  such  a  large  fre¬ 
quency  range  that  intrinsic  structure  might  dominate.  If  scattering 
is  important,  we  expect  [3  ss  —2.  We  find  an  average  value  of  /3  = 
—  1.9  ±  0.1  for  the  entire  sample.  Clearly,  we  anticipate  that  in¬ 
trinsic  structure  may  be  important  for  some  sources,  but  that  the 
(3  is  close  to  the  expected  value  from  scattering  indicates  that 
scattering  is  important  for  the  sample  of  sources  as  a  whole. 

Table  1  summarizes  the  inferred  scattering  and  intrinsic  diam¬ 
eters  from  fitting  the  data  for  each  source  to  equation  (1).  For 


comparison,  we  also  list  the  predicted  diameter  at  1  GHz  from 
the  NE2001  model  (Cordes  &  Lazio  2002)  for  interstellar  scatter¬ 
ing,  0NE2OO1  ■  Figure  1  illustrates  examples  of  the  measured  angular 
diameters,  showing  the  results  for  an  AGN  for  which  a  relatively 
large  amount  of  scattering  is  inferred  and  one  for  which  a  rela¬ 
tively  small  amount  of  scattering  is  inferred. 

There  are  values  of  the  reduced  Table  1  that  are  larger 
than  unity,  at  times  by  a  significant  factor.  These  result  from  a 
combination  of  two  factors.  First,  the  uncertainties  for  some  data 
are  likely  to  be  underestimated.  For  the  observed  diameters  ob¬ 
tained  from  our  observations,  we  estimated  their  uncertainties 
using  a  bootstrap-like  procedure  in  which  visibility  data  asso¬ 
ciated  with  different  antennas  were  removed  before  performing 
the  fit.  The  range  of  fitted  diameters  suggests  a  1 0%  uncertainty. 
For  angular  diameters  obtained  from  the  literature,  we  have  as¬ 
sumed  the  same  fractional  uncertainty.  However,  we  have  been 
able  to  identify  data  for  which  this  assumed  1 0%  is  likely  to  be 
too  small. 

The  fitting  procedure  also  yields  an  estimate  of  the  uncertainty 
in  the  inferred  scattering  diameter  (at  1  GHz).  The  median  value 
of  this  uncertainty  is  0.1  mas.  We  have  repeated  the  fitting  pro¬ 
cedure  with  larger  uncertainties  adopted  for  the  angular  diameters 
obtained  from  the  literature,  and  in  some  cases  even  removing  ap¬ 
parent  outliers.  The  typical  change  in  the  inferred  diameter  is  com¬ 
parable  to  the  uncertainty  in  the  inferred  scattering  diameter.  We 
have  also  repeated  analyses  described  below  with  larger  uncer¬ 
tainties  for  the  angular  diameters  from  the  literature  and  find  no 
change  from  the  results  we  present  below. 

A  second  potential  cause  of  large  x  ^  values  is  that  we  have 
adopted  fixed  frequency  scaling  exponents  in  equation  (1).  We 
might  obtain  a  better  fit  by  letting  x  be  a  fitted  parameter,  fitting 
the  scattering  frequency  dependence  rather  than  adopting  —2.2, 
or  both.  To  do  so  would  often  require  more  data  than  are  avail¬ 
able.  Consequently,  while  larger  than  unity,  we  consider  there  to 
be  plausible  explanations  for  the  x^  values  in  Table  1  and  will  use 
the  angular  diameters  resulting  from  our  fits. 
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TABLE  2 

Statistical  Measures 

Mean/Median  Scintillators  Nonscintillators 


Galactic  Latitude  (Absolute  Value) 


Mean  (deg) .  29.5  ±  3.0  46.3  ±  7.0 

Ecliptic  Latitude  (Absolute  Value) 

Mean  (deg) .  34.9  ±  4.3  28.8  ±  4.9 

Scattering  Diameters  (at  1  GHz) 

Mean  (mas) .  2.1  ±  0.2  2.1  ±  0.2 

Median  (mas) .  1.9  2.0 

Redshifts 

Mean  (mas) .  1.38  ±  0.18  1.29  ±  0.31 

Median  (mas) .  1.27  1.40 


Ojha  et  al.  (2006)  divided  the  sources  into  two  groups,  “scin¬ 
tillators”  and  “nonscintillators.”  Since  the  publication  of  that 
paper,  it  has  been  realized  that  some  of  the  AGNs  identified  as 
scintillators  displayed  variation  at  only  a  single  epoch  (“once-er”) 
among  the  MASIV  observations,  leading  to  the  possibility  that  a 
nonscintillator  would  be  classified  mistakenly  as  a  scintillator. 
Analysis  of  the  MASIV  observations  (to  be  published  elsewhere) 
suggests  that  the  light  curve  from  an  individual  epoch  can  be  clas¬ 
sified  correctly  at  the  95%  confidence  level.  From  the  four  epochs 
of  observations  comprising  MASIV,  the  probability  of  a  false  iden¬ 
tification  is  only  4.3%,  meaning  that  we  expect  only  one  source 
in  our  sample  to  be  classified  mistakenly.  While  we  identify  the 
single-epoch  variable  sources  in  Table  1 ,  their  presence  should 
have  no  significant  efifect  on  our  analysis,  and  we  treat  the  single¬ 
epoch  variable  AGNs  as  scintillators.  Also,  analysis  has  continued 
on  the  MASIV  sources,  so  there  may  be  occasional  differences  in 
the  classification  (scintillating  vs.  nonscintillating)  in  our  Table  1 , 
as  compared  to  Table  1  of  Ojha  et  al.  (2006). 

Figure  2  shows  the  distribution  of  the  inferred  scattering  diam¬ 
eters  plotted  as  a  function  of  Galactic  coordinates.  We  have  con¬ 
sidered  the  distribution  of  the  sources  on  the  sky  as  a  function  of 
both  Galactic  latitude  b  and  ecliptic  latitude  (3  and  the  distribution 
as  a  function  of  solar  elongation.  There  is  no  statistically  signifi¬ 
cant  correlation  of  the  inferred  scattering  diameter  with  either  co¬ 
ordinate  (strictly,  the  absolute  value  of  the  coordinate),  nor  with 
cos  (|h|)  or  cos  (|/3|).  (Typical  correlation  coefficients  are  approx¬ 
imately  0.1.)  The  use  of  the  cosine  of  the  angle  [cos(|h|)  or 
cos  (|/3|)]  attempts  to  compensate  for  the  increased  amount  of 
sky  at  low  latitudes.  There  is  also  no  correlation  between  inferred 
scattering  diameter  and  solar  elongation. 

The  lack  of  a  correlation  of  the  scattering  diameter  with  Galac¬ 
tic  latitude  would  appear  to  be  at  odds  with  previous  determina¬ 
tions  that  interstellar  scattering  increases  rapidly  at  low  latitudes 
(e.g.,  Rao  &  Ananthakrishnan  1984).  However,  while  no  longi¬ 
tude  selection  criterion  was  applied  in  constructing  our  sample, 
our  sample  is  nonetheless  weighted  strongly  toward  the  outer 
Galaxy.  We  have  no  sources  with  longitudes  in  the  range  —60°  < 
/  <  60°,  and  only  a  few  at  low  latitudes  in  the  range  —120°  < 
/  <  120°.  Thus,  we  attribute  the  apparent  lack  of  a  correlation 
between  scattering  diameter  and  Galactic  latitude  as  a  result  of 
having  few,  essentially  no,  lines  of  sight  into  the  inner  Galaxy. 

While  there  is  no  correlation  of  scattering  diameter  with  Ga¬ 
lactic  latitude  for  the  entire  set  of  sources,  scintillating  AGNs  are 


Scattering  Diameter  (mas) 

Fig.  3. — Histograms  of  scattering  diameters.  The  open  histogram  shows  the 
distribution  for  the  scintillating  sources;  the  hatched  histogram  shows  the  dis¬ 
tribution  for  the  nonscintillating  sources. 

consistently  at  lower  Galactic  latitudes  than  the  nonscintillating 
AGNs.  Table  2  shows  that  the  average  absolute  values  for  the  Ga¬ 
lactic  latitudes  of  scintillating  and  nonscintillating  sources  differ 
by  nearly  20°.  No  such  difference  is  found  for  the  average  (abso¬ 
lute)  ecliptic  latitude. 

Figure  3  shows  a  histogram  of  the  inferred  scattering  diameters. 
We  have  used  a  Kolmogorov- Smirnov  (K-S)  test  to  assess  whether 
scintillating  AGNs  have  a  different  distribution  of  inferred  scatter¬ 
ing  diameters,  as  compared  to  the  nonscintillating  AGNs.  We  find 
no  statistical  difference:  the  scattering  diameters  of  scintillating 
AGNs  do  not  differ  appreciably  from  those  of  the  nonscintillating 
AGNs. 

Examination  of  the  scattering  diameters  inferred  for  individ¬ 
ual  objects  indicates  that  some  of  the  largest  scattering  diameters 
result  from  AGNs  for  which  no  measurements  exist  below  1  GHz. 
In  order  that  these  not  bias  our  result,  we  removed  these  and  re¬ 
peated  the  K-S  test  analysis.  There  is  no  change  in  the  result,  that 
the  scattering  diameters  for  the  scintillating  and  nonscintillating 
AGNs  are  consistent  with  having  been  drawn  from  the  same  dis¬ 
tribution.  Both  the  mean  and  the  median  scattering  diameter  for 
scintillating  and  nonscintillating  sources  is  approximately  2  mas 
(Table  2). 
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Redshift 

Fig.  4. — Distribution  of  scattering  diameters  as  a  function  of  redshift.  Stars 
indicate  sources  that  scintillate;  squares  indicate  nonscintillating  sources. 
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From  the  entire  sample,  37  AGNs  (64%)  have  measured  red- 
shifts.  There  appears  to  be  little  difference  in  the  redshift  distri¬ 
bution  of  the  scintillating  and  nonscintillating  AGNs,  with  the 
two  populations  having  similar  means  and  medians  (Table  2). 
Figure  4  shows  the  distribution  of  the  scattering  diameters  as  a 
function  of  redshift.  We  have  determined  the  correlation  between 
the  scattering  diameters  and  redshifts  for  the  entire  sample,  as  well 
as  splitting  it  into  the  two  populations,  scintillating  and  nonscin¬ 
tillating.  There  is  no  correlation  of  the  scattering  diameter  with 
redshift  for  the  entire  sample.  There  may  be  a  marginal  correla¬ 
tion,  at  the  5%  confidence  level,  between  the  scattering  diameters 
and  redshift,  in  the  opposite  sense  for  the  scintillating  and  non¬ 
scintillating  sources.  That  is,  the  scattering  diameters  of  scintil¬ 
lating  (nonscintillating)  AGNs  may  become  smaller  (larger)  at 
higher  redshifts. 

4.  DISCUSSION  AND  CONCLUSIONS 

In  one  sense,  our  results  are  broadly  consistent  with  what  is 
known  about  intraday  variability  (IDV)  and  interstellar  scatter¬ 
ing.  In  our  sample,  scintillating  AGNs  lie  typically  at  lower  Ga¬ 
lactic  latitudes,  consistent  with  the  notion  that  the  scintillation 
responsible  for  IDV  results  from  the  Galactic  ISM  (§  1). 

In  addition  to  the  ISM,  other  ionized  media  along  the  lines  of 
sight  to  these  sources  are  the  interplanetary  medium  (IPM)  and 
the  intergalactic  medium  (IGM).  It  is  possible  that  the  IPM  could 
contribute  angular  broadening  at  a  level  sufficient  to  be  detectable 
in  this  analysis,  particularly  given  the  use  of  0.33  and  0.61  GHz 
observations.  The  lack  of  a  correlation  with  either  the  ecliptic  lat¬ 
itude  (Table  2)  or  solar  elongation,  however,  implies  that  the  IPM 
makes  no  detectable  contribution  in  these  observations.  The  lack 
of  a  strong  correlation  with  redshift  and  the  systematically  lower 
Galactic  latitude  of  the  broadened  sources  suggests  that  the  IGM 
makes  no  detectable  contribution  to  our  observations,  either. 

In  a  further  effort  to  differentiate  the  Galactic  contribution  of 
scattering  from  any  possible  intergalactic  contribution,  we  have 
searched  for  pulsars  within  1°  of  the  AGNs  in  our  sample.  We  find 
no  pulsars  this  close  to  any  of  our  sources.  Given  the  relatively  low 
density  of  pulsars  on  the  sky,  a  significantly  larger  sample  of  AGN  s 
would  be  required  in  order  to  make  such  a  comparison. 

On  the  face  of  it,  the  result  that  scintillating  AGNs  are  broad¬ 
ened  at  levels  comparable  to  those  of  the  nonscintillating  AGNs 
appears  to  contradict  the  requirement  (§  1 )  that  in  order  to  display 
scintillation,  a  source  must  be  sufficiently  compact.  One  possi¬ 
bility  is  that  scattering  is,  in  fact,  not  important  at  all,  and  that  we 
have  mistakenly  attributed  the  effects  of  intrinsic  structure  to  an¬ 
gular  broadening.  We  regard  this  as  unlikely,  given  that  the  an¬ 
gular  diameters  of  the  sample  of  sources,  as  a  whole,  scale 
approximately  as  expected  from  interstellar  scattering. 

In  fact,  the  magnitude  of  the  estimated  broadening  does  not 
appear  to  be  problematic  from  the  standpoint  of  quenching  the 
scintillations.  Taking  2  mas  at  1  GHz  as  a  characteristic  scatter¬ 
ing  diameter  (Table  2),  the  implied  scattering  diameter  at  5  GHz 
is  80  /ias,  comparable  to  the  value  derived  for  a  number  of  the 
well-known  scintillating  sources  (e.g.,  Rickett  et  al.  1995,  2002; 
Dennett-Thorpe  &  de  Bruyn  2003;  Bignall  et  al.  2003). 

A  notable  feature  of  the  most  extreme  IDV  sources  is  that  de¬ 
tailed  analyses  of  their  light  curves  suggest  that  the  scattering 
medium  responsible  for  the  scintillation  lies  quite  close  to  the 
Sun  (ss25  pc).  If  local  material  was  responsible  for  the  scintil¬ 
lation  of  all  scintillating  AGNs,  we  would  not  expect  a  difference 
between  the  average  Galactic  latitude  of  scintillating  and  non¬ 
scintillating  sources  (Table  2).  That  such  a  difference  exists  sug¬ 
gests  that  the  scintillation  for  most  scintillating  AGNs  results  from 
scattering  material  associated  with  the  Galactic  ISM. 


A  consistent  explanation  for  these  results  is  obtained  if  the  scin¬ 
tillation  is  produced  from  small  “clumps”  of  scattering  material, 
distributed  throughout  the  Galactic  disk.  For  instance,  Dennett- 
Thorpe  &  de  Bmyn  (2003)  determine  that  a  “thin  screen”  is  re¬ 
sponsible  for  the  extreme  IDV  of  J1 8 19-1-3845,  with  the  screen 
being  about  1 0  pc  distant  and  having  an  internal  level  of  scatter¬ 
ing  measured  by  =  0.5  They  do  not  provide  a  quan¬ 

titative  estimate  of  the  thickness  of  this  screen,  but,  in  order  that 
the  screen  be  “thin,”  it  must  surely  be  the  case  that  its  thickness 
is  AL  <  1  pc. 

Suppose  we  require  that  the  scattering  contributed  by  such  a 
clump  not  make  a  significant  contribution  to  the  angular  broad¬ 
ening.  For  illustration  purposes,  we  adopt  the  amount  of  broad¬ 
ening  contributed  by  the  clump  to  be  ci  "^0.2  mas,  which  would 
be  only  a  1 0%  contribution  to  the  typical  broadening  diameter  that 
we  measure.  The  resulting  scattering  measure  (SM;  Cordes  & 
Lazio  2002)  is  then  SMd  ~  10^"*  ^  kpc  If  the  clump  has 

cl  =  0.5  the  implied  thickness  is  0.05  pc  (~10"^  AU), 

which  would  certainly  qualify  as  thin.  Moreover,  following  Rickett 
et  al.  ( 1 995 )  and  Rickett  (2002),  it  can  be  shown  that  a  more  distant 
scattering  clump  tends  to  produce  a  lower  scintillation  index. 

Thus,  our  picture  is  one  in  which  the  Galactic  disk  contains 
(small)  clumps  of  scattering  material.  Lines  of  sight  through  the 
disk  are  scattered  by  the  broadly  distributed  ionized  ISM,  so  that 
AGNs  over  the  range  of  latitudes  that  we  observe  have  measur¬ 
able  scattering  diameters.  Some  (many?)  lines  of  sight  pass 
through  one  of  these  clumps,  and  AGNs  having  compact  enough 
components  are  then  observed  to  scintillate.  However,  the  clumps 
are  small  enough  that  they  produce  effectively  no  additional  broad¬ 
ening.  This  scenario  is  also  broadly  consistent  with  the  notion  of 
clumps  of  material  producing  extreme  scattering  events  (Fiedler 
et  al.  1987)  and  parabolic  arcs  in  pulsar  dynamic  spectra  (Hill 
et  al.  2005). 

We  can  also  use  the  difference  between  the  scintillating  and 
nonscintillating  sources  to  set  quantitative  limits  on  the  amount 
of  radio  wave  scattering  contributed  by  the  IGM.  We  adopt 
0.5  mas  at  1  GHz  («3  a  from  Table  2)  as  the  upper  limit  on  the 
difference  in  the  amount  of  scattering  between  the  two  popula¬ 
tions.  The  implied  scattering  measure  is  SM  ^  10^^  kpc 
(Cordes  &  Lazio  2002).  In  turn,  the  SM  is  given  by 

SM  =  Csm7^£>,  (2) 

where  D  is  the  distance,  F’  is  a  fluctuation  parameter  encapsulat¬ 
ing  aspects  of  the  microphysics  of  the  plasma,  rig  is  the  electron 
density,  and  Csm  =  1.8  cm®  is  a  constant.  For  a  character¬ 

istic  redshift  of  approximately  unity  (Fig.  4),  the  equivalent  (angu¬ 
lar  size)  distance  is  D  «  1.5  Gpc,  implying  Fnl  %  10^'®  ®  cm^®. 

For  a  baryonic  matter  density  =  0.127  (Spergel  et  al. 
2007),  the  intergalactic  electron  density  can  be  no  larger  than 
We  <  2.2  X  10^^  cm^^,  assuming  that  helium  is  fully  ionized 
(Sokasian  et  al.  2002).  Thus,  we  require  F  ^  10^,  so  as  not  to 
violate  the  inferred  limits  on  scattering.  For  reference,  in  the  dif¬ 
fuse  Galactic  ISM,  F  «  0.2,  and  in  the  Galactic  spiral  arms, 
F  ^  10.  In  turn,  the  F’-parameter  is 

F  =  CeS-Hfl\  (3) 

where  is  the  normalized  second  moment  of  the  fluctuations,  e 
is  the  fractional  variance  in  «g  within  the  plasma,  rj  is  the  filling 
factor,  and  /o  is  the  largest  scale  on  which  the  density  fluctua¬ 
tions  occur  (or  outer  scale,  if  the  plasma  is  turbulent),  in  parsec 
units.  Assuming  that  C  ^  e  ~  1,  we  conclude  that  ^  10^^. 
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The  IGM  is  thought  to  be  permeated  by  shocks  (Dave  et  al. 
200 1 ),  which  might  be  expected  to  drive  77  ^  1 .  Given  the  larger 
scales  available  in  the  IGM,  /q  1  Mpc  would  not  be  unreason¬ 
able.  We  are  forced  to  conclude  that  the  current  limits  on  inter- 
galactic  scattering,  while  broadly  consistent  with  the  current 
picture  of  the  IGM,  do  not  yet  place  significant  constraints  on  its 
properties. 

While  we  find  no  indications  of  intergalactic  scattering,  future 
observations  are  warranted.  In  particular,  if  a  scintillating  AGN 
is  found  close  to  the  line  of  sight  to  a  pulsar,  a  comparison  be¬ 
tween  the  two  lines  of  sight  would  provide  strong  constraints  on 
the  amount  of  Galactic  versus  intergalactic  scattering.  Also,  higher 
sensitivity  observations  (e.g.,  with  the  VLBA  High  Sensitivity 
Array  [HSA])  targeting  scintillating  AGNs  with  larger  diameters 
may  provide  additional  constraints.  Many  of  the  AGNs  with  the 
largest  diameters  are  not  detected  at  the  lower  frequencies,  frequen¬ 
cies  at  which  the  VLBA  alone  has  a  relatively  low  sensitivity.  The 
HSA  could  be  used  to  verify  whether  these  AGNs  do  indeed  have 
such  large  scattering  diameters,  or  assess  to  what  extent  intrinsic 
structoe  contaminates  the  scattering  diameter  estimates. 

We  summarize  our  findings  as  follows.  In  our  sample  of 
58  AGNs,  approximately  75%  of  the  sample  exhibit  IDV  (inter¬ 
stellar  scintillation),  with  the  other  25%  not  showing  IDV.  Inters¬ 
tellar  scattering  is  measurable  for  most  of  these  AGNs,  and  the 
typical  broadening  diameter  is  2  mas.  Scintillating  AGNs  are  typ¬ 
ically  at  lower  Galactic  latitudes  than  the  nonscintillating  AGNs, 
consistent  with  the  scenario  that  IDV  is  a  propagation  effect  from 
the  Galactic  ISM.  The  magnitude  of  the  inferred  interstellar  broad¬ 
ening  measured  toward  the  scintillating  AGNs,  when  scaled  to 
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higher  frequencies,  is  comparable  to  that  determined  from  anal¬ 
yses  of  the  light  curves  for  the  more  well-known  intraday  variable 
sources.  However,  we  find  no  difference  in  the  amount  of  scatter¬ 
ing  measured  toward  the  scintillating  versus  nonscintillating  AGNs. 
A  consistent  picture  is  one  in  which  the  scintillation  results  from 
localized  regions  (clumps)  distributed  throughout  the  Galactic 
disk,  but  that  individually  make  little  contribution  to  the  angular 
broadening.  In  our  sample,  63%  of  the  AGNs  have  measured 
redshifts.  At  best,  a  marginal  trend  is  found  for  scintillating  (non¬ 
scintillating)  AGNs  to  have  smaller  (larger)  angular  diameters  at 
higher  redshifts.  Finally,  while  broadly  consistent  with  the  sce¬ 
nario  of  a  highly  turbulent  IGM,  our  observations  do  not  place 
significant  constraints  on  its  properties. 
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